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ABSTRACT 

We present a catalog of Faraday rotation measures (RMs) and redshifts for 4003 extragalactic radio sources 
detected at 1 .4 GHz, derived by identifying optical counterparts and spectroscopic redshifts for linearly polar- 
ized radio sources from the NRAO VLA Sky Survey. This catalog is more than an order of magnitude larger 
than any previous sample of RM vs. redshift, and covers the redshift range < z < 5.3; the median redshift of 
the catalog is z = 0.70, and there are more than 1500 sources at redshifts z > 1. For 3650 of these sources at 
Galactic latitudes \b\ > 20°, we present a second catalog in which we have corrected for the foreground Faraday 
rotation of the Milky Way, resulting in an estimate of the residual rotation measure (RRM) that aims to isolate 
the contribution from extragalactic magnetic fields. We find no significant evolution of RRM with redshift, 
but observe a strong anti-correlation between RRM and fractional polarization, p, that we argue is the result 
of beam depolarization from small-scale fluctuations in the foreground magnetic field or electron density. We 
suggest that the observed variance in RRM and the anti-correlation of RRM with p both require a population 
of magnetized intervening objects that lie outside the Milky Way but in the foreground to the emitting sources. 
Keywords: catalogs — galaxies: distances and redshifts, magnetic fields — magnetic fields — polarization 



1. INTRODUCTION 

Faraday rotation is a powerful probe of magnetic fields 
along the line of sight between a linearly polarized radio 
source and the observer. When combined with redshift infor- 
mation, Faraday rotation measurements provide the potential 
to constrain the evolution of cosmic magnetic fields over time 
scales corresponding to the redshifts of the most distant polar- 
ized sources detectable. However, previous studies of Faraday 
rotation measure (RM) vs. redshift have been hampered by the 
small size of the available samples (200-300 sources at most), 
often drawn from highly inhomogeneous data sets. Radio and 
optical surve ys such as the NRAO Very Large Array Sky Sur- 
vey (NVSS; ICondon et aflll998l) . and the Sloan Digital Sky 
Survey (SDSS: lYork et al. 2000) can enable a comprehensive 
new analysis of RM vs. redshift, thereby allowing much more 
sensitive studies of the evolution of magnetic fields over cos- 
mic time. 

1.1. Faraday Rotation 

Faraday rotation is a birefringence effect whereby the linear 
polarization angle {8) of a radio wave is rotated when prop- 
agating through a magnetized, ionized gas. For a distant lin- 
early polarized radio source observed at a wavelength A, the 
polarization angle is rotated by an amount A(9 = RM A 2 , where 
RM is defined as: 



RM(z. s ) = 0.81 



' n e (z)B\\(z) dl 
(l+z) 2 dz 



dz rad m 
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In this Equation, z s is the redshift of the polarized source; n e (z) 
is the free electron number density at some foreground red- 
shift z, in cm" 3 ; B||(z) is the line-of-sight component of the 
magnetic field at redshift z, measured in fiG; and dl is a line- 
element along the line of sight, measured in parsecs. A pos- 
itive RM corresponds to a magnetic field oriented toward the 



observer. In the idealized case of a Faraday thin system, one 
can determine the RM by measuring values for 8 at multiple 
observing wavelengths, and then performing a linear fit to 8 
vs. A 2 . 

Equation (Q]i implies that the observed RM is a super- 
position of contributions from multiple magneto-ionic re- 
gions along the entire line of sight. In every case, there is 
a contribution from the Ear th's ionosphere (^1-2 rad m" 2 ; 
iBrentjens & de Bruynl l2005h . and for extragalactic sources 
there is always a foreground contribution from the Milky 
Way known as the Galactic Rotation Measure (GRM; ^10- 
1000 rad m~ 2 , varying with Galactic longitude and latitude; 
ISchnitzelerl 120 1 Ot lOppermann et alj 120121) . Outside of the 
Galaxy, there are many contributions to the observed RM, in- 
cluding the polarized sources themselves, intervening galax- 
ies, their halos, intracluster gas or independent clouds of mag- 
netized gas. The extragalactic component of RM can be stud- 
ied by calculating the residual rotation measure (RRM), such 
that RRM = RM - GRM. 

1.2. Previous Catalogs of RM and Redshift 

Over the past several decades, there have been numerous 
studies in which RRMs of polarized sources as a function 
of redshift have been us ed to probe extragalactic magnetic 
fields over cosmic time dBernet et alJ 1200a : iFuiimoto et al ' 
Kronberg & Simard-Normandin 1976; Kronberg etal 
Kronberg & Perrvlll982t iKronberg et al.l l2008h iNelson 



& Wolfd 119951; iReinhardtl 119721; ISofueetal. 
;on & Nelsonfl 982: Wel ter et alJll984tlYou et ~ 



1973: lOren 
1979; Thomson 
2003). These studies have demonstrated that the mean RRM 
of sources over the sky is zero out to redshifts z > 4. How- 
ever, the variance of the RRM distribution, the effect of in- 
tervenors/absorbers and the precise nature of any RRM evo- 
lution with redshift all remain contentious. The largest sam- 
ple reported to date is the (as-yet unpublished) data set ana- 
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lyzed by Kron berg et alj d2008f) . whose sample consisted of 
268 RRM measurements for sources out to z ~ 3.7. These 
authors presented evidence for a significant increase in the 
variance of the RRM distribution as a function of increasing 
redshift, from which they proposed that galaxies hosted strong 
magnetic fields at relatively early times. 

In the present paper, we use archival surveys and databases 
to derive a new catalog of RM and redshift for more than 4000 
extragalactic radio sources out to redshifts z ~ 5.3, a sample 
that is markedly larger and extends to much higher redshifts 
than all previous efforts. We also present an additional catalog 
in which we have corrected for Galactic Faraday rotation for 
most sources, resulting in a compilation of RRM vs. z suitable 
for studies of the evolution of Faraday rotation and magnetic 
fields as a function of cosmic time. In ^2] we describe the 
databases and surveys used in our analysis. In |3] we present 
our catalog of RM vs. z, while in $4] we present a catalog of 
RRM vs. z- In J5]we discuss the overall properties of the sam- 
ple, and consider the relationships between RRM, fractional 
polarization and redshift for these sources. 

2. DATABASES AND SURVEYS 

In the following subsections, we describe the radio and op- 
tical data sets from which we obtained values of RM and z, 
respectively. The sky coverages of the surveys used in this 
study are indicated in FigureQ] the final catalog contains only 
sources with declinations S > -40°, constrained by the cover- 
age of the NVSS (from which all our RMs have been derived). 

2.1. Obtaining Rotation Measures 

The NVSS dCondon etalJI 19981) was a 1.4 GHz radio sur- 
vey that used the Very Large Array (VLA) to image the entire 
sky for <5 > -40° at an angular resolution R w 45". The result- 
ing source catalog contains information on Stokes parameters 
/, Q and U for over 1.8 million discrete objects. 

The RMs that we use for our analysi s come from the re - 
processing of the original NVSS data bvlTavlor et alj d2009l) . 
and are available onlineQ Tayl or et all d2009l) used images of 
Stokes Q and U in two frequency bands at 1364.9 MHz and 
1435. 1 MHz for every NVSS source to determine linearly po- 
larized intensity P = \JQ 2 + U 2 and fractional linear polariza- 
tion p = P/I. If the signal-to-noise ratio for P was greater than 
8, they derived the RM from the difference in polarization po- 
sition angle across the two bands. The resultant distribution 
of 37 543 RMs is provided in Figure [2] Because the RMs 
were derived only from two closely spaced frequency bands, 
the data are potentially subject to nn ambiguities and non-A 2 
effects. However, subsequent b roadband polarization studies 
for subsets of the iTaylor et aD d2009t) catalog have produced 
res ults that are in up to 96% agreement with the RMs derived 
bvlTavlor et al.l J2009h dMao et al.ll2010t iVan Eck et al.ll20TTt 
lLaw et al.ll201ll) . 

2.2. Obtaining Redshifts 

Our sources of redshift data are two online databases plus 
several recent large optical redshift surveys. We exclude from 
consideration any objects flagged as stellar or Galactic. We 
also only use spectroscopic redshifts (as opposed to photomet- 
ric redshifts or Lyman dropouts) due to the greater reliability 
and sufficient abundance of spectroscopic redshifts from the 
surveys we cons ider. D etails o n these databases and surveys 
are provided in {32.2.1l through $2.2.61 

^http : //www . ucalgary . ca/ras/rmcatalogue| 



To assign redshifts to the RM catalog of Taylor et al. 
diool . we need to make associations between sources de- 
tected at both optical and radio wav elengths. There have been 
numerous approac h es to this issue (jMagliocchetti et al.l l 1998 ; 
l lvezic etal.1 12002[ iMcMahon ~t all 120021: iBest et all 2005 : 
Sadler et all 20071: iKimball & Ivezid 120081: iLin etafl boiO : 
Plotki net al.l |2Q10T) . Such studies use angular proximity as 
the primary discriminant in making an association between 
an optical source and a radio source, although the proximity 
requirement varies widely from 2" to 3', depending on the 
priority placed by different investigators on finding all gen- 
uine associations or minimizing false matches. Some authors 
have considered the possibility of multi-component sources 
with varying degrees of complexity, from merely incorporat- 
ing double-lobed radio sources to consid eration of all of dou- 
ble, t r iple, bent and co re-jet scenarios (Maglio cchetti et aTl 
119981: IBest et al.l 120051) . Matching procedures also differ 
based on whether associations are made with a relatively low- 
resolution radio survey such as the NVS S (R = 45"), a high- 
resolution survey such as FIRSlQ (R = 5": lBecker etaL| | 1995 ; 
Ivezic et al.l2002l). or a combination of both (Best et al] |2005 : 
Plotki net al.l l2010f) . Some procedures make use of proper- 
ties s uch as source size, flux or inclination be tween radio 
lobes (Mag liocchetti et al.ll998LlBest et al.l2005l) . while other 
methods calculate exc ess matches above those from random 
sourc e distributions (Kimball & Ivezic 120081: IPlotkin et al.l 
l2010h . 

2.2.1. NED 

The NASA/IPAC E xtragalactic Database (NED; 
iHelou et all 119911 I1995IFI stores data on >170 million 
objects outside the Milky Way. NED includes 4 million 
sources with redshift information, including many associ- 
ations between radio sources and optical redshifts. These 
data come from surveys, published results and references to 
known redshifts in the literature. We discarded any redshift 
that was not obtained spectroscopically, and also discard 
entries for which the "z Quality" field was flagged. 

2.2.2. SIMBAD 

The Set of Identifications, Measure ments and Bibliogra - 
phy for Astronomical Data (SIMBAD; IWenger et all [2000 fl 
is a database operated out of the Center de Donnees as- 
tronomiques de Strasbourg (CDS) in Strasbourg, France. It 
contains nearly five million Galactic and extragalactic objects. 
Like NED, SIMBAD includes redshift data for radio objects 
derived from surveys and from the literature. 

2.2.3. SDSS 

The SDSS3 dYork et al.ll2000h is an optical survey that has 
so far imaged more than one-third of the sky (14500 square 
degrees) in five different wavebands at R ks 0.5", using a ded- 
icated 2.5 m optical telescope at Apache Point, New Mex- 
ico. For the analysis in this paper w e use data release 8 
(DR8), released on 2011 January 11 (lAihara et al.l 1201 ll) . 
which contains spectroscopic data (including redshifts) for 

2 The V LA Faint Images of the Radio Sky at Twenty centimeters survey 
(see j|2.3l for more details). 

3 |http : //ned . ipac ■ caltech . edu/| , accessed on 2012 August 16. 

4 |http : //simbad . u- strasbg. f r/ simbad/| accessed on 
201 2 August 17. 

5 |http: / /sky server . sdss3 ■ org/dr8/en/| 
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approximately 860000 galaxies and 116000 quasars over al- 
most 9300 square degrees. We used Structured Query Lan- 
guage (SQL) to obtain data on spectroscopic redshifts and 
optical magnitudes for galaxies and quasars from the online 
DR8 database. We applied filters to our search to select only 
objects classified as galaxies or quasars, and to selectively 
remove any sources which had warnings associated with the 
redshift measurements. 

2.2.4. The 6dFGS Survey 

The Six-d egree Field Galaxy Survey (6dFGS; iJones et ail 
2004, 20090 measured the spectroscopic redshifts of 1 10 256 
galaxies over 17000 square degrees, with a median redshift 
Z = 0.053. The survey used the 6dF fibre-fed multi-object 
spectrograph on the UK Schmidt telescope of the Australian 
Astronomical Observatory (AAO). We applied filters to our 
search to select only objects in the top two redshift quality 
bands, as determined by the survey (quality 3 or 4). 

2.2.5. The IdFGRS Survey 

The Two-degree Field Galaxy Redshift Survey (2dFGRS; 
IColless et all 1200 lL 120031(1 was carried out with the 2dF in- 
strument on the Anglo-Australian Telescope (AAT). It cov- 
ered approximately 1500 square degrees and obtained reliable 
redshifts for 221 414 galaxies at a median redshift z = 0.11. 
We applied filters to select only objects in the superior red- 
shift quality bands, as determined by the survey (quality 3-5). 

2.2.6. The 2QZ/6QZ Survey 

The 2dF QSO Redshi ft survey (2QZ) an d the 6dF QSO 
Redshift survey (6QZJ3 ( Croom et al.| [2004) were the quasar 
spectroscopic redshift survey counterparts, respectively, to 
the 2dFGRS and 6dFGS described above. In combination 
they are referred to as the 2QZ/6QZ survey. Combined, the 
2QZ/6QZ survey contains redshift data for 23 660 quasars 
over 720 square degrees. We applied filters to select only ob- 
jects with high-quality redshifts, as determined by the survey 
team (Quality_Flag_l = 11, 21 or 31). 

2.3. Morphological Information 

We additionally make use of the higher angular resolution 
of the FIRST survey, with which we match radio sources 
with optical d ata using morphol ogical information. The 
FIRST survey (Becker et al. 1995) consists of VLA images 
at 1.4 GHz covering 9900 deg 2 , matched to the sky coverage 
of SDSS. FIRST detected approximately one million sources, 
and its high angular resolution (R = 5" compared to R = 45" 
for the NVSS) enables the detection of compact sources and 
provides detailed morphological information. 

3. THE RM-REDSHIFT CATALOG 

Previous catalogs o f RM v s. z such as those pre- 
sented bv IWelter et all (fl98l . lOren & Wolfel (fl995l) and 
iKronberg et al.1 (120081) have suffered from a high degree of in- 
homogeneity. Typically, both RM and z data have been com- 
piled from a range of small published (and occasionally un- 
published) samples, which were then collated into an RM-z 
catalog. This approach was undertaken to maximize the size 
of the data set, although the results were still only catalogs 

6 http : //www-wf au ■ roe .ac.uk/ 6dFGS/| 

7 http : //www2 . aao ■ gov ■ au/2dFGRS/| 

8 http : //www. 2df quasar . org/Spec_Cat/ c atalogue . html 



of 100-200 sources. With the wealth of new data available, 
we here attempt to strike a balance between creating a cata- 
log with a large number of sources, while also maintaining a 
high degree of homogeneity in the data sets from which our 
inf ormation comes. T o associate optical redshift data with 
the lTaylor et alJ (120091) RM catalog, we first identify sources 
for which associations had previously been made in NED or 
SIMBAD, an d the n use data from the optical redshift surveys 
described in 32.2l to make new associations. 

3.1. NED and its Limitations 

NED is currently the largest single database of extragalactic 
sources, and therefore represents the natural starting point to 
search for redshifts. A total of 2023 of the redshifts that we 
ultimately include in our RM-redshift catalog were provided 
from information in NED (50.5% of 4003 RM-redshift pairs). 
These represent cases for which a specific NVSS source has 
an associated redshift within the NED data base. 

Though many of the Taylor et al. ( 2009) NVSS sources are 
listed within NED, testing with the Sesame Name Resolvefl 
revealed that some NVSS sources are missing. Correspon- 
dence with the NED team (M. Schmitz, 2011, private com- 
munication) confirmed that only around 90% of the NVSS 
source catalog is in NED. The database is complete for right 
ascensions h < RA < 5 h 40 m , but outside this range the cov- 
erage is incomplete. Inquiries revealed that source inclusion 
is being deferred until associations are made with preexisting 
entries in NED. This task has been ongoing for three years, 
with no firm completion date planned. Those sources not 
included have complications such as a single NVSS source 
being resolved into several FIRST sources or the presence 
of more than one object in the relevant field. Unfortunately, 
those sources not included are precisely those that are espe- 
cially relevant for this project, as they are radio galaxies that 
have likely also been observed in the optical by spectroscopic 
surveys. We ther efore make these associations ourselves, as 
described in 33.21 

3.2. A Radio-Optical Association Algorithm 

We add to the redshifts for NVSS sources provided by 
NED and S I MBAD by associating the radio sources in the 
iTaylor etafl (120091) catal og w ith redshifts from the four opti- 
cal surveys described in 32.21 There are, however, a number 
of difficulties with making associations between optical and 
radio data. One significant example is galaxies which have 
complex radio morphologies, such as double-lobed structures, 
that are not mimicked in the optical. To help make these as- 
sociations, we include data from the higher resolution FIRST 
survey (though some NVSS sources either lack a FIRST de- 
tection or lie outside the FIRST sky coverage). There is 
good complementarity between the high angular resolution 
of FIRST and the superio r sensitivity of NVSS to extended 
emission dBest e t al. 2005), allowing detection of both com- 
pact and multi-component radio sources. 

Our association a lgorithm is primarily based on the work 
of iBest et alJ J2005), who tested and refined their criteria us- 
ing Monte Carlo simulations involving catalogs of random 
sky locations. The association criteria arising from these 
simulations, based on the entire NVSS catalog of 1.8 mil- 
lion sources, were preferred to attempting to formulate new 
criteria using simulations on the much smaller subset of 

9 |http : //cdsweb . u-strasbg . f r/cgi-bin/Sesame| 



4 



Hammond, Robishaw & Gaensler 



37 543 source listed by Tayl or etaT] d2009|). We also incor - 
porate aspects of t he sc hemes of Magliocch etti et all dl998l) . 
Ilvezic et all d2002l) . and iMcMahon et alj d2002l) ; a summary 
of the main elements of our algorithm is represented in Fig- 
ure[3] Excluding associations taken from NED and SIMBAD, 
the algorithm contributes 89% of the RM-z pairs in our RM- 
redshift catalog; the remainder were flagged by the algorithm 
as having complex features and were subsequently manually 
inspected. 

3.2.1. Association Classes 

Our algorithm associates sources and assigns them to one 
of five association classes, A, B, C, D or E, as summarized 
in Table [T] For sources in class A, we have used FIRST data 
to make the association; for sources in class B, we rely on 
NVSS dat a onl y (the particular limitations of class B are de- 
scribed in 33.5b ; class C corresponds to more complex sources 
that required manual visual inspection; class D corresponds to 
associations made within NED, while class E indicates asso- 
ciations made by SIMBAD. Class A is further subdivided into 
seven subclasses, A(i) through A(vii), depending on the radio 
morphology of the source in the FIRST data. 

The admission criteria for association classes A to C are fur- 
ther explained below. Examples of sources in each of classes 
A(i), A(ii), A(iii), A(iv), A(v), A(vi), A(vii), B and C are pro- 
vided in Figuresg]|3|6l|7][|]|il[I21[II]and[T2] respectively. 

Following iBest et all (|2P05), we account for the possibil- 
ity of multi- component NVSS sources by making an initial 
search of the lTaylor et alJ (2009) catalog in a 3' radius around 
each optical source. The distance of 3' is selected to en- 
sure that genuine multi-component sources have at least two 
matches, while still being much less than the typical separa- 
tion of 8'- 10' between NVSS sources. The chance of random 
matches within 3' is high, motivating more complex selection 
criteria as described below. As shown in Figure[3] our scheme 
automates associations in cases where there are fewer than 
three NVSS matches. 

In cases for which there is one NVSS source within 3' of 
a given optical source, our algorithm checks for FIRST data 
within 30" of the optical position. If FIRST data are found 
and if the NVSS-optical offset is less than 15", the algorithm 
proceeds. A single NVSS source more than 15" from the op- 
tical source is not considered an association. The next steps 
taken by the algorithm depend on how many FIRST matches 
are found within 30". 

A single FIRST match is accepted if the FIRST-optical off- 
set is less than 3" [class A(i); see examples in Fig. |4j| or if the 
FIRST-optical match is within 10" and is also within 75% of 
the major axis of the optical source [class A(ii); Fig.|3. The 
rate of false detections for matches within 3" is estimated to 
be less than 1%, while the inclusion of matches within 10", 
subject to the additional criterion for class A(ii), helps main- 
tain high completeness. 

A double or triple match within FIRST is accepted if one 
of the two radio sources is within 3" of the optical galaxy 
[classes A(iii) and A(v), shown in Figs.|6]and[8] respectively] 
or if the doubles criterion is met for two of the matched FIRST 
sources [classes A(iv) and A(vi), shown in Figs. [7] and [9] re- 
spectively). These require that the quantity O x S x F is less 
than five, where O is the arcsecond offset between the optical 
galaxy and the mean position of two FIRST sources, S is the 
ratio of source size for those two FIRST sourc es, and F is their 
flux ra tio. This derives from the finding of McMahon et al. 
(2002), that double-lobed radio sources detected in FIRST 



have two lobes of similar sizes and fluxes, and a mean po- 
sition that is close to the optical galaxy. The reliability of this 
procedure is estimated to be > 99%. 

Our scheme also treats the case where two NVSS sources 
are found within 3' of the same optical source. These sources 
are classified as potential double-lobed radio systems if a set 
of offset criteria are met: both radio sources must be within 
90" of the optical galaxy, the mean position of the two ra- 
dio sources must match the optical position within 15", and 
the nearer of the two radio components to the optical source 
must be > 15" from the optical position (or there is less than 
20" between the two radi o components). W e then search for 
FIRST data within 30". Bes t et all d2005l) did not consider 
sources in the absence of FIRST data, but noted the likely de- 
gree of contamination in the case of do ubles if NVSS data 
alone are relied upon. For this work, the iTaylor et"aT] (12009b 
catalog contains sufficiently few candidate doubles that those 
without clear FIRST data (i.e., a FIRST source within 3" of 
the optical source position) are visually inspected (see be- 
low). Potential doubles with clear FIRST data are associated 
as doubles and placed in class A(vii). Note that sources in 
class A(vii) have two NVSS sources (the two lobes of a ra- 
dio galaxy) matched to the same optical source. Within the 
RM-redshift catalog, the two components of each pair are ar- 
bitrarily designated class A(vii)a or A(vii)b in reflection of 
this. Figure [10] shows two examples of class A(vii) associa- 
tions. 

In some cases for which there is only one NVSS match, 
there are no corresponding FIRST data because the source 
lies outside the sky coverage of the FIRST survey (see Fig- 
ure Q]), or because the source has faded between the NVSS 
and FIRST observations, or because a radio source seen as 
substantially extended in NVSS is resolved out in the higher 
angular resolution FIRST data. In these cases, a match is ac- 
cepted if the NVSS-optical offset is less than 10", and the 
source is then placed in class B. Figure [TT] shows two ex- 
amples of class B sources. The choice of a 10" matching 
radius is conservative in light of other work, and avoids sig- 
nificant contamination from background so urces (B est et al] 
l2005HPlotkin et alJ2010l:ISadler et al.l2007l) . It is also smaller 
than the NVSS-optical offset limit of 15" that we use for class 
A (i.e., when a FIRST counterpart is detected). 

For sources for which there are three or more NVSS 
matches, four or more FIRST matches, or a candidate dou- 
ble source with no clear FIRST data, the algorithm flags the 
optical source as requiring manual visual inspection. In these 
complex situations, images of the relevant parts of the sky 
are downloaded from the NVSS and FIRST postage stamp 
servers, the SDSS image server and the SuperCOSMOS im- 
age extraction service (for 2dFGRS, 6dFGS and 2QZ/6QZ). 
In cases for which we confirm visually that the optical and ra- 
dio sources are associated, then t he value for redshift from the 
optical survey is assigned to the ITaylor etaT] (120091) source, 
and these sources are assigned class C. Two examples of class 
C sources are shown in Figure [12] While visual inspection 
introduces some level of subjectivity to the resulting associ- 
ations, we note that only 195 associations in our final RM- 
redshift catalog fall into this class, just 4.9% of the total. 

3.2.2. Noteworthy Features of the Association Algorithm 

In 25 cases (23 of which were drawn from 6dFGS), the al- 
gorithm matches two different sources f rom the same optica l 
survey with a single radio source from iTaylor eTail (2009). 
In all such cases, we visually inspected the optical and radio 
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images, using morphology and proximity to determine which 
optical source is associated with the radio source. 

We note that while our association scheme attempts to use 
NVSS and FIRST together to achieve the best possible com- 
promise between completeness and reliability, our approach 
does exhibit a slight bias against extended sources. This is be- 
cause nearly all matches involving a single-component FIRST 
source are detected by our algorithm, while not all of those in- 
volvin g multiple NVSS components are identified (B est et all 
l200l . There is no a priori reason that this slight bias should 
impact the results of our RM analysis. 

3.3. Choosing a Final Redshift for the RM-Redshift Catalog 

Many of the polarized sources in the iTavlor etail d2009l) 
catalog have an optical counterpart and resulting redshift in 
more than one of the optical surveys or databases that we 
have considered. We list every such redshift in our RM- 
redshift ca talog, but additiona lly identify a "selected redshift" 
for each ITavlor et all (120091) source, representing that red- 
shift which we consider to be most reliable. For example, 
59 NVSS sources have redshifts from two separate optical 
surveys, while one source was associate d with three optica l 
surveys. In 2296 cases, associations of ITavlor et al.l (2009) 
sources are made with bo th databases and surveys (for ex- 
ample, a ITavlor et al.l 120091 source that it is associated with 
an SDSS optical counterpart by our algorithm will sometimes 
have the same NVSS/SDSS match listed in NED), or an op- 
tical counterpart is listed both in NED and in SIMBAD. In 
many cases, the decision on a selected redshift is trivial, ei- 
ther because all redshift entries are the same (730 cases), or 
because all the associated redshifts agree closely with one an- 
other (1439 cases with < Az < 0.01). From our total of 
4003 sources, this leaves 127 cases for which multiple dis- 
tinct redshifts are found, so that the selection of a final red- 
shift could have a functional impact upon the catalog. The 
differing redshifts are due to disagreements between the con- 
stitutive surveys/databases. Some of these disagreements are 
attributable to identifiable errors such as measuring the red- 
shift of an absorber along the line of sight to an emitter and 
wrongly assigning that redshift to the emitter. 

In those cases where a substantive redshift selection needs 
to be made, the final selected redshift is chosen from among 
detections in multiple optical surveys/online catalogs in the 
descending priority listed in Table [2] This prioritization gives 
preference to surveys over databases, and prefers surveys with 
high angular resolution and wide-area sky coverage. NED is 
preferred to SIMBAD because of the greater number of reli- 
able redshifts in NED. 

3.4. Contents of the RM-Redshift Catalog 

As described in ^2.21 and ^3] we have obtained r edshifts 
for the polarized radio sources of Taylor et al. (2009) from a 
range of optical databases and surveys, with each optical/radio 
match given a corresponding association class as defined in 
Table Q] The resulting RM-redshift catalog contains 4003 
sources with both RMs and redshifts (with a "selected" red- 
shift listed in cases where multiple values for z have been 
found). This catalog contains an order-of-magnitude more 
entries than any previous compilation of RM vs. redshift, and 
extending out to much higher redshifts (z = 5 .27). A sample of 
rows and columns from the RM-redshift catalog is shown in 
Table |3j the full RM-redshift catalog (available in its entirety 
as an online machine-readable table) contains the coordinates, 



fractio nal polarization, R M and all other data for each source 
from lTavlor et al.l ( 120091) . along with positions, redshifts, pho- 
tometry, object types, journal references and association class 
derived for each optical match, and the selected redshift and 
other relevant data associated with the best-matching optical 
counterpart. 

Table]4]shows the contents of the RM-redshift catalog bro- 
ken down by survey and association class. The most common 
association class is A(i), representing a simple unresolved ra- 
dio source that aligns closely with the corresponding optical 
counterpart. This is followed by class B, representing a simi- 
lar situation to class A(i) but in cases either where no FIRST 
source is detected or no FIRST data are available. The 2dF- 
GRS, 6dFGS and 2QZ/6QZ surveys provide most of their as- 
sociations in class B, as their sky coverage has little overlap 
overlap with FIRST. In contrast, almost all SDSS associations 
fall in class A, because the sky-coverages of FIRST and SDSS 
were intentionally matched. The total number of associations 
arising from each survey varies as expected according to the 
sky coverage and number of objects in each survey. For ex- 
ample, most survey associations are made with the SDSS, 
covering 9300 square degrees and containing almost 1 mil- 
lion extragalactic redshifts, while the fewest associations are 
made with 2QZ/6QZ, containing only 23 000 quasar redshifts 
and focusing on two small regions of the sky. The row la- 
be led "Total" in Table |4]indicates the total number of matches 
to ITavlor et"afl (120091) for each survey or database, while the 
column labeled "All" indicates the total number of times each 
association class is assigned to a radio/optical match. The 
rows and columns labeled "Selected" show the contribution 
from each survey/database and from each association class, 
respectively, to the final selected redshifts using the priority 
order listed in Table [2] 

3.5. Completeness, Reliability and Testing of Matching 
Processes 

It is desirable to estimate both the completeness and the re- 
liability (also known as efficiency; cf. Kimball & Ivezic l2008h 
of our RM-redshift catalog. Completeness refers to the frac- 
tion of real matches found, while reliability refers to the frac- 
tion of matches that are real. Both these quantities are difficult 
to estimate for our RM-redshift catalog because they are re- 
liant upon a number of different factors: the completeness and 
proportion of spurious sources within each constituent survey, 
the number of true matches expected, and the level of back- 
ground contamination in different regions of the sky. 

Previous studies involving matches of NVSS with FIRST 
have found high completen e ss and reliability dBest et alJ 
l200HlKimball & Ivezic12008l) . iKimball & Ivezicl (120081) con- 
sidered a range of matching radii; for the radius of 30" that we 
adopt, they found a completeness of 99.7% and a reliability of 
96%. Matching between FIRST and SDSS yields a complete- 
ness above 94% and reliability above 90%. This rises further 
when NVSS data are included, with iBest et alJ (120051) sug- 
gesting 95% as a conservative completeness estimate. 

Figure [13] shows the difference in NVSS and optical sky 
positions for matched radio-optical sources. The distributions 
of offsets in both Right Ascension and Declination are sharply 
peaked around zero, with a small spread. More than 90% of 
matched NVSS sources lie within 3" of the associated optical 
source; when double-lobed radio sources are excluded, this 
fraction rises to 94%. 

One particular issue is the 375 sources in the RM-redshift 
catalog in association class B, i.e., sources that were asso- 
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ciated by our algorithm using optical and NVSS data only, 
without a FIRST detection. Class B associations are likely 
to be less re liable than those that also draw on FIRST data. 
Specifically, iBest et afl ((2005) have used Monte Carlo simu- 
lations to estimate that this association class may have a false 
detection rate of up to 6%. However, class B represents less 
than 10% of the entries in the RM-redshift catalog. Remov- 
ing these sources would cause the completeness of the catalog 
to suffer, and also would further bias the sample against ex- 
tended sources (due to the removal of extended sources that 
lie within the FIRST sky coverage but lack a FIRST detec- 
ti qn). For the RM-r edshift catalog as a whole, the simulations 
of lBest et all d2005l) suggest that schemes that use both NVSS 
and FIRST have an overall reliability > 98%. 

We tested our algorithm using the p rior radio-optical asso- 
ciation work of | Oren & Wolfd d 19951) and iKimball & Ivezid 
(2008). lOren & Wolfd (11995b provided a small test sample 
of 20 sources that we used to test the basic functionality of 
o ur scheme. A more rob ust test involves the larger catalog 
of IKimball & Ivezid (120081) . who associated the SDSS with a 
range of rad i o surv eys, including the NVSSH Of the 37 354 
iTavlor et alJj 2009) sources, 1 295 appear with redshift data in 
the wo rk of (Kimball & Ivezid (120081) . 

The IKimball & Ivezid d2008l) catalog used SDSS DR6 
( Adelman-McCarthv et al. 2008), so we also performed tests 
using that data set, downloaded fr om the SDSS archive. Us - 
ing the matching radii selected by Kimball & Ivezid (|2008), 
which were different from our own due to the differing aims 
of that study, we were able to mak e 126 of the 1295 asso- 
ciations made by Kimball & Ivezic (2008). The small differ- 
ences are attributable t o different correlation programs (the 
Kimball & Ivezic 2008 scheme, for instance, takes the FIRST 
data, not the optical data, for the position of the source), 
rounding differences arising from c omputation using d i fferen t 
programs, and oddities within the Kimball & Ivezid (2008) 
catalog. 

We also tested the robustness of our approach by randomly 
selecting 100 sources that had been associated via our algo- 
rithm, and visually inspecting them. In all cases, the visual 
inspection confirmed the association made by the code. This 
indicates that our algorithm associates optical and radio de- 
tections with a reliability > 99%. 

4. THE RRM-REDSHIFT CATALOG 

4.1. Subtraction of the Galactic RM Contribution 

In considering the evolution of magnetic fields over cosmic 
time, it is insufficient to simply plot RM against redshift be- 
cause the foreground Faraday rotation from our own Milky 
Way must first be accounted for. This is especially important 
because the Galactic contribution is not uniform, but varies 
across the sky. 

Due to the increased spatial fluctuations in the GRM at low 
Galactic latitudes, the first step in calculating RRM is to ex- 
clude all sources wi th \b\ < 20° (draw ing upon the study of 
these fluctuations by Schnitzeler 2010, and taking a threshold 
for \b \ co mparable to that chos en for studies by Welt er et"aT] 
119841 and lOren & Wolfdll995l) . We thus only calculate an 
RRM for 3650 of the 4003 sources presented in CT4l the 
breakdown of which for each optical survey/database and for 
each association class is listed in the final row and final col- 
umn of Table |U respectively. 

10 The Kimball & Ivezic 12008) catalog is available online at 

http : / / www .astro. Washington, edu/user s/akimball/racliocat 



A range of methods have been employed in the literature 
for calculating and subtracting the GRM. In this paper, we 
use the new map of the f oreground Faraday sky computed by 
Opperman n et al.l (120121) . as shown in Figure [14] This map 
enables the calculation of an RRM by subtracting the GRM at 
that sk y posit ion from the RM va lue reported by Tayl or et al.l 
(2009). The ITavlor eTafl (12009ft RM data were used in the 
construction of the lOpperman n et al.l (120121) map, but the re- 
construction algorithm used in that work filters out the ex- 
tragalactic contributions, ensuring that the resul ting map is 
still ap propriate to use in GRM calculations for Tayl or et al.1 
(2009) sources. 

Figure Q3] compares RM and RRM for our data. The 
efficacy of this GRM subtraction is indicated by the con- 
siderable narrowing of the RRM distribution (mean +0.3 ± 
0.4 rad m~ 2 , standard deviation 23.2 rad m" 2 ) compared to 
the RM distribution (mean +2.3 ±0.6 rad m~ 2 , standard de- 
viation 36.1 rad m" 2 ). We also considere d a number of 
previously published RRM schemes (lOren & Wolfd 119951: 
Johnst on-Hollitt et all l2004t ISofue et all 119791: IShortetalJ 
120071) . but the GRM map of Opper mann et al.l d2012l) resulted 
by far in the smallest scatter in RRM of all approaches con- 
sidered. 

4.2. Contents of the RRM-Redshift Catalog 

After applying the GRM correction discussed in jii4.ll we 
derive a catalog of 3650 radio sources with both RRMs and 
redshift. An extract of some columns from this RRM-redshift 
catalog can be seen in Table[5] The full RRM-redshift catalog 
(available online as a machine-readable table), contains the 
coordinates, flux, fractional polarization, RM, GRM, RRM, 
redsh ift, plus appropriate ancillary data for the radio source 
from ITavlor etalj d2009) and for the best-matching optical 
counterpart. 

This RRM-redshift catalog is more than an order of mag- 
nitude larger than any such catalog previously published. 
This is highlighted by Figure [16] in which we compare 
our catalo g data against the samples of RRM vs. z de- 
scribed bvlkronberg & Perrvl dl982l) . [Welter et al.l dl98l and 
lOren& Wolfd (119951) . 

5. ANALYSIS OF THE RRM-REDSHIFT CATALOG 
5.1. Summary of Source Characteristics 

The RRM-redshift catalog of 3650 sources covers a redshift 
range from to 5.27, with a median redshift z = 0.70 and in- 
cluding almost 1400 sources with redshifts z > 1. The RRM- 
redshift catalog covers an RM range from -465.4 to +270.4 
rad m~ 2 , and an RRM range from -476.5 to +206.1 rad m~ 2 . 
The RRM distribution, however, is tightly clustered around 
zero: 88% of sources have |RRM| < 25 rad m" 2 , 96% have 
I RRM < 50 rad m" 2 , and 99.3% have |RRM| < 100 rad irT 2 . 
The 1 .4 GHz Stokes / fluxes of the sample ranges from 1 1 mJy 
to 55 Jy, with a median of 300 mJy. Polarized fluxes range 
from 2.6 mJy to 1.3 Jy, with a median of 8 mJy. 

The RRM-redshift catalog contains a diverse range of 
sources, drawn from several surveys and databases, and with 
a range of multi-wavelength properties and environments. A 
full analysis of this data set is beyond the scope of this pa- 
per, but here we present a brief analysis of its overall prop- 
erties. To address potential issues of inhomogeneity in the 
catalog, we consider the full data set plus three relatively 
large and well-defined subsets: the 1376 sources for which 
/ihe selected redshift was drawn from SDSS (referred to as 
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"SDSS sources" in subsequent discussion), the 516 sources 
with SDSS redshifts for which the optical spectrum was clas- 
sified by SDSS as a galaxy ("SDSS galaxies"), and the 860 
sources with SDSS redshifts for which the optical spectrum 
was classified by SDSS as a quasar ("SDSS quasars"). 

For each of these four data sets, we consider three relation- 
ships: RRM as a function of z; fractional linear polarization, 
p, as a function of z; and RRM as a function of p. We note 
that any relation between RRM and p does not explicitly re- 
quire knowledge of the redshift of the sources as has been the 
focus of this paper. However, because we have optical iden- 
tifications, consideration can be given to how the RRM vs. p 
relation changes for different subsets of the overall data such 
as SDSS galaxies vs. SDSS quasars. 

5.2. Residual Rotation Measure vs. Redshift 

Figure [T7] shows the distribution of RRM as a function of 
redshift for (a) the entire RRM-redshift catalog (this is the 
same data as shown in the bottom panels of Figs. Q3] and 
[TBI (b) all SDSS sources, (c) SDSS galaxies, and (d) SDSS 
quasars. As expected, there is a marked difference in the red- 
shift distribution of SDSS galaxies (median redshift z = 0.17) 
compared to SDSS quasars (median redshift z = 1 .22). 

As indicated by the red lines overlaid on each panel of Fig- 
ure [17] there is no apparent trend in the mean or standard de- 
viation of RRM as a function of z in any of the four sam- 
ples consideredPl This is in strong contrast to the results of 
IRronberg et al.1 (120081) . who presented a significant increase 
in the standard de v iation of RRM to higher redshifts. 

IRron berg et al] (12008b further characterized an evolution 
of RRM with z by splitting their data into two groups, cor- 
responding to sources above and below a threshold redshift 
Zb- By applying a Kolmogorov-Smirnov (KS) test to the nor- 
malized cumulative distri butions of |RRM| at redshifts above 
and below Zb, IRronberg et alj ( 120081) found that the RRMs of 
low- and high-redshift sources differed at 99% significance 
for zi, ~ 1-8. We can repeat this experiment with far larger 
sample sizes: in our RRM-redshift catalog, there are 3140 
and 510 sources at Zb < 1 -8 and Zb > 1 -8, respectively. Apply- 
ing a KS test to the distributions of |RRM| on either side of 
this threshold, the values of |RRM| for the two samples differ 
at 53% significance, which is consistent with the two groups 
of data being drawn from the same underlying distribution. 

As a more sensitive test of any evolution of RRM with z, 
we use the Spearman rank test to look for evolution of |RRM| 
as a function of z for all four data sets shown in Figure[T7] We 
find no correlation for panel (c), and find a weak correlation 
between |RRM| and z, at 2.5er, 2.2er and 3.8cr significance for 
panels (a), (b) and (d), respectively. We do not consider any 
of these trends significant, especially given that here and in 
further subsections below we are examining a diverse range 
of different possible correlations. 

We thus conclude that there is no significant dependence of 
RRM or its variance as a function of z in our RRM-redshift 
catalog, in contrast to the strong effect o f this kind reported 
previously using much smaller data sets dWelter et al.l 1 1984c 
lYou et alj|2003b IRronberg et al.ll2008l) . Recently. iBernet etalj 

1 1 The two places in Figure [T7l a) where the standard deviation is substan- 
tially larger than ^20 — 25 rad m are produced entirely from two sources 
with anomalously large RRMs (a source at z = 1 .44 with RRM = -444 rad irT 2 
and a source at z = 2.13 with RRM = —476 rad itT 2 , respectively). The latter 
source is an SDSS quasar, and so also has an effect on the standard deviations 
in panels (b) and (d). 



d2012|) have considered a much smaller sample of 371 RMs 
from iTavlor et a l. (2009), and reported a sim ilar inability to 
reprod uce the RRM vs. z behaviour seen by Kronberg et al. 
(2008). Given this lack of redshift dependence, we now 
briefly consider possible origins of the observed RRMs. 

We first note two terms that must contribute to the stan- 
dard deviation of 23.2 rad m" 2 seen for the RRMs in Fig- 
ure [TTfa): the measu rement errors of the individual RMs in 
the lTavlor et al.l (120091) catalog, and the error asso ciated with 
the GRM calculation of lOpperman n et ail (120121) . The stan- 
dard deviation in RM due to measur ement errors is 1 1 rad m" 2 
(ISchnitzeleri 120 lOt IStil et alJ 1201 ll) . while the mean error 
in GRM for the sources in our RRM-redshift catalog is 
6 rad m" 2 . Subtracting these in quadrature from the observed 
variance leaves a standard deviation of 20 rad m" 2 that must 
come from one or more astrophysical phenomena. 

An obvious explanation for the lack of redshift evolution 
is that the RRMs are solely a residual contribution from the 
Milky Way, resulting from imperfect foreground subtraction, 
or from RM fluctuations on small er angular scales than ar e 
being sampled by the GRM map of Op permann et alJ (120121) . 
We can rule out the former option, because while the RMs 
of our sources show a strong dependence on Galactic longi- 
tude and latitude (as seen in Fig. [2] for the larger NVSS RM 
sample from which our catalog is derived), the corresponding 
RRMs show no pattern or trend with Galactic coordinates. To 
consider the latter opt ion, we note that the typ ical spacing be- 
tween the sources that Opper mann et al.l d2012l) used to calcu- 
late their foreground map is ~ 1°, meaning that any Galactic 
contribution to our observed RR Ms must represe nt GRM fluc- 
tuations on angular scales < l°. IStil et alJ(l201 ll) calculate the 
fluctuations in GRM on these scales and show that, at latitudes 
\b\ > 20° as we are considering here, the standard deviation in 
RM on a scale ofl°is 12-17 rad m" 2 . Subtracting this Galac- 
tic term in quadrature from our observed variance implies that 
there is a 10-15 rad m" 2 contribution to the standard devia- 
tion in RRM that cannot be due to small-scale fluctuations in 
GRM, and hence must be extragalac tic. This can be compared 
with the work of Schnitzeler (2010), who performed a st atis- 
tical decomposition of the RMs in the lTavlor et alJ d2009l) cat- 
alog into different components as a function of Galactic lat- 
itude, and concluded that extragalactic Faraday rotation con- 
tributed a standard deviation of « 6 rad m" 2 to the RMs of 
TayloretaD (120091) . An agree ment within a facto r of ~ 2 be- 
tween our estimate and that of lSchnitzelerl (l2010h seems rea- 
sonable given the very different approaches taken between the 
two studies. 

We now consider the possibility that most of the extragalac- 
tic contribution to RRM arises in the polarized sources them- 
selves, e.g., in an envelope of ionized gas in the host galaxy 
or in its immediate environment. If the RRM originates at the 
same redshift as the emitter, then we expect a (1 +z)~ 2 dilution 
factor as per Equation ([JJ). The corresponding contribution 
to the variance in RRM should then decrea se with redshift 
in Fig ure [171 (cf. the dashed line in Fig. 7 of Kronber g et alJ 
2008). Specifically, if we assume that the extragalactic con- 
tribution to the standard deviation in RRM of 10-15 rad m" 2 
corresponds to sources at our median redshift z = 0.7, then an 
identical population of sources at a redshift z = 2 should only 
contribute a standard deviation in RRM of 1-2 rad m" 2 . Con- 
sidered in quadrature with the other terms contributing to the 
scatter in RRM as discussed above, this dilution should result 
in an overall decrease in the standard deviation of RRM by 



8 



Hammond, Robishaw & Gaensler 



2-4 rad m~ 2 between z = 0.7 and z = 0.2. Although this is 
a small effect, we can rule out its presence in our data — 
the Spearman rank test discussed above shows that if any- 
thing, the variance in RRM slightly increases, not decreases, 
with redshift. The lack of evolution of RRM with z can then 
only be explained if the standard deviation in RRM evolves 
with (1 +z) 2 in the emitter's reference frame, to cancel out the 
( 1 +z)~ 2 effect as the signal propagates to Earth. This would be 
a fortuitous coincidence and would also represent very strong 
evolution of the emitted RRM with redshift, both of which 
make this possibility unlikely. 

The alternative is that the extragalactic component of RRM 
towards each source is introduced in one or more intervening 
systems along the line of sight between the polarized source 
and the observer, a poss ibility considered in detail by many 
previous authors (e . g.. iBernet et al.l 120081: I Kronberg et alj 
l2008t IQren & Wolfdll995t IWelter et aUll9841) . If the inter- 
vening systems are all at comparable redshifts, then we will 
not see a (1 +z)~ 2 dilution term, and indeed might expect a 
slight increase in variance of RRM with increasing z, since 
more distant sources are more likely to have an intervenor 
along the line of sigh t (see solid and dotted lines in Fig. 7 of 
IKronberg et~ai1l2008l) . IKronberg etaT] (120081) advocate a pop- 
ulation of intervenors with a standard deviation of RRM of 
~60-l 15 rad m" 2 in the observer's frame, but our analysis ar- 
gues that a contribution 10-15 rad m" 2 is more likely. This is 
broadly consistent with the halo of a Milky- Way-like galaxy 
(n e » 3 x 10" 4 cm" 3 , B ps 1 pG, I » 50 kp c; iGaensler et all 
20081 iMao et afll2"010l ISun & Reichl l2012). but a proper in- 
terpretation requires further, detailed consideration of this and 
other possible intervening source populations. 

5.3. Polarized Fraction vs. Redshift 

Figure [18] shows how p varies with z for the entire RRM- 
redshift catalog, for all SDSS sources, for SDSS galaxies, and 
for SDSS quasars. Panels (a) and (b) show a clear trend, in 
that sources with z < 0.5 can have a range of polarized frac- 
tions extending beyond p = 30%, while sources at z > 0.5 are 
confined solely to lower polarized fractions, p ~ 5%. 

Panels (c) and (d) provide a simple explanation for this ap- 
parent bimodality. Within the SDSS sub-sample, the low- 
redshift sources that can be both weakly or strongly polarized 
are all classified as galaxies based on their optical spectra. In 
contrast, the SDSS sources that extend to high redshifts and 
that are all relatively weakly polarized all have optical spectra 
indicating that they are quasars. While this distinction needs 
further study, a likely explanation is that the NVSS counter- 
parts to the SDSS galaxies are radio lobes from active galax- 
ies, which are expected to be Faraday thin and thus show high 
degrees of polarization. In contrast, the NVSS counterparts 
to SDSS quasars are radio-loud cores, for which optical depth 
effects and strong magnetic fields are expected to result in re- 
duced polarization levels. 

We can confirm this hypothesis by considering the distribu- 
tion of association classes for the SDSS galaxies compared to 
the SDSS galaxies. As summarized in Table Q] the most ro- 
bust classifications are class A, in which an NVSS source has 
a clear match with an optical counterpart and also shows a 
relatively simple morphology in FIRST. Within class A, sub- 
classes A(i) and A(ii) are characteristic of core-dominated ra- 
dio morphologies, while classes A(iii) through A(vii) repre- 
sent radio morphologies for which a significant fraction of 
the emission is from radio lobes. If the radio polarization from 



SDSS galaxies is mainly from radio lobes and that from SDSS 
quasars is mainly from radio cores, then we expect the SDSS 
galaxies to be dominated by classes A(iii) to A(vii), while 
the SDSS quasars should be dominated by classes A(i) and 
A(ii). Indeed the data support this expectation: considering 
only SDSS sources in class A, 35% of SDSS galaxies are in 
classes A(i) and A(ii), while 65% are in classes A(iii) through 
A(vii). In contrast, 54% of SDSS quasars are in classes A(i) 
and A(ii), while 46% are in classes A(iii) to A(vii). Fur- 
thermore, we note that class C corresponds to complex, ex- 
tended morphologies that cannot be classified automatically. 
For SDSS galaxies, 20% of all sources in the RRM-redshift 
catalog are in class C, while for SDSS quasars, the fraction 
of class C sources drops to just 8%. All this additional mor- 
phological information clearly demonstrates that SDSS galax- 
ies are dominated by lobes, while SDSS quasars are domi- 
nated by cores, consistent with the differing distributions of 
fractional polarization for these two populations seen in Fig- 
ureUll 

Noting the differing polarization properties of galaxies and 
quasars, we can then consider whether there is any evolution 
of p with z in Figures \Wic) and (d). Applying the Spearman 
rank test, we find no significant correlation of p with z for 
either SDSS galaxies or SDSS quasars, as apparent by eye via 
the red lines marking the mean and standard deviation of p as 
a function of z. A similar absen ce of evolution of p w ith z at 
1 .4 GHz was recently found by Banfiel d'et al.l d201 lh from a 
smaller sample of 69 sources in the redshift range 0.04 < z < 
3.2. 

Other studies of polarization as a function of redshift 
have usually considered parameters such as the "depolariza- 
tion measure" (DP; the ratio of p at a high frequency to 
that at a low frequency) or X^n (the wavelength at which 
p falls to half its peak value) (e.g., IKronberg et al.l 11972b 
iGoodlet et alJl2004l) . Many of these studies have found that 
DP increases with z or that A1/2 decreases with z, in both 
cases indicating increased depolarization at higher redshifts 
dConwav et al]| 19741: IGoodlet et all l2004t IGoodlet & Kaiser! 
l2005Hlnoue & Tabaral 19821: IKronberg et alj|1972l) . While we 
do not see such an effect here, it is important to note that such 
studies convert observed values of DP or A1/2 to those seen 
in the reference frame of the emitter. We cannot make such 
a correction for our sample because we do not have multi- 
wavelength information that would tell us how p depends on 
w avelength for e a ch sou rce. With the exception of the study 
of iBanfield eTai] (f20TTh . our result is thus not directly com- 
parable with previous work on the depolarization of radio 
sources as a function of redshift, without the inclusion of sev- 
eral additional assumptions. 

5.4. Residual Rotation Measure vs. Polarized Fraction 

In Figure [19] we plot the dependence of RRM on p for the 
entire RRM-redshift catalog, for all SDSS sources, for SDSS 
galaxies, and for SDSS quasars. In this case we see a striking 
pattern in all four panels: the standard deviation in RRM is 
large for low fractional polarizations, and is small for high 
fractional polarizations. 

We first note that this is not an artifact of sensitivity and 
signal-to-noise. The uncertainty in an individual RM mea- 
surement (and hence the variance in a large sample of similar 
RM measur ements) is inversely proporti onal to the signal-to- 
noise ratio (Brentiens & de Bruvn 2005). Thus we expect a 
source with a low polarized flux to have a large error in RM, 
and vice versa. However, we are here plotting the fractional 
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polarization rather than the polarized flux, and there is no cor- 
relation between these two parameters in our catalog. (If any- 
thing, any bias is in the reverse direction, since the faintest 
sources can only be detected in polarization if their fractional 
polarization is high.) The behavior seen in Figure |T9]is thus 
astrophysical, rather than instrumental. 

The most likely interpretation for this behavior is a de- 
polarization mechanism, since stronger Faraday effects in 
or in front of a polarized source can induce re duced polar- 
ization through a variety o f mechanisms (e . g.. iBurnl Il966t 
Gardn er & Whiteoakl Il966t iTribbld 11991b ISokoloff et all 
1998). Bandwidth depolarization cannot be a contributor 
in Figure [19] since this effect is only significant for RMs 
with magnitudes l arger than ~ 100 rad m" 2 (see Fig. 1 of 
iTavlor et al.l 120091) . Just 2% of the sources in our RRM- 
redshift catalog have |RM| > 100 rad m" 2 , so this cannot be a 
significant contributor to the observed depolarization. Depth 
depolarization is also unlikely to be a factor, since this only 
occurs when the emitting medium is mixed with t he m edium 
producing the Faraday rotation. We discussed in §15.21 above 
how the lack of evolution of RRM with z argues against an 
intrinsic origin for the RRMs, meaning that the RRM signal 
occurs wholly in the foreground to the polarized sources and 
so cannot produce depth depolarization. 

The remaining possibility is beam depolarization, whereby 
high values of RRM imply large fluctuations in RRM on 
scales smaller than the angular extent of the radio source 
(even in cases where the source is unresolved by the telescope 
beam). Such fluctuations will cause differing polarization an- 
gles within the beam to cancel, resulting in a reduced polar- 
iz ed fraction. 

lHaverkorn et al.l (|2008) have shown that the Galactic fore- 
ground can indeed simultaneously produce larger RMs and 
enhanced beam depolarization against unresolved polarized 
background sources. Beam depolarization due to the Galactic 
foreground produces an anti-correlation between GRM and p 
(lHaverkorn et al.l 120081) . and so thus can possibly produce a 
dependence between RRM and p given that one component 
of RRM is likely to be a residual Galactic cont ribut ion pro- 
duced by small-scale fluctuations in GRM (see 35.21 i. While 
the depolarizing effects of the Galaxy are likely to be present 
in our data, they are insufficient on their own to explain the 
trend that we see bet ween |RRM| and polariz ed fraction. For 
example, Figure 6 of lHaverkorn et alT f2008) shows that for 
similar angular resolution and observing frequency to those of 
the NVSS data being considered here, fluctuations in Galac- 
tic RM at the level of ^200 rad m" 2 are needed to depolar- 
ize background sources by a factor of ~ 2. In contrast, Fig- 
ure [19] suggests that fluctuations in RRM only of magnitude 
~ 50 rad m" 2 can depolarize background sources by a fac- 
tor of ~3-5. The depolarization that we observe is thus too 
strong an effect to be explained only by small-scale Galactic 
RM fluctuations. The beam depolarization also does not seem 
to be originating in th e emitting sources themselves, since we 
have argued in ^5.21 above that the RRMs do not have a sig- 
nificant intrinsic component. A significant contribution to the 
observed depolarization must therefore be due to small-scale 
fluctuations in RM somewhere between the source and the 
Milky Way, in intervening magneto-ionic material along the 
line of s ight. This is the same origin independently proposed 
in ^5.21 for the extra g alacti c component of the RRMs them- 
selves. iBernet et al.l (120121) have used a different approach, 
focusing on the frequency dependence of depolarization, but 



have separately come to a similar conclusion that the extra- 
galactic intervenors contributing to the observed RRMs must 
be highly turbulent. 

6. CONCLUSIONS 

We h ave compare d the e xtragalactic rotation measure cat- 
alog of ITavlor et aLl (120091) with redshift data from a range 
of optical surveys and databases to produce a new sample 
of 4003 sources with both RM and redshift data (the "RM- 
redshift catalog"). We further derive a catalog of residual 
rotation measure vs. redshift, in which we have subtracted 
the Galactic RM contribution toward a subset of 3650 high- 
latitude sources (the "RRM-redshift catalog"). The result- 
ing samples contain more than an order of magnitude more 
sources than any previously published catalogs of RM vs. red- 
shift or RRM vs. redshift. 

In our RRM-redshift catalog, we do not see any signifi- 
cant evidence that the variance of RR M changes wi t h red- 
shif t, in contrast to previo us studies by IWelter et al.l (119841) 
and lKronberg et all Q008-) who found that the RRMs of then- 
sources showed increased scatter at higher z. The overall stan- 
dard deviation of the RRMs in our catalog is 23 rad m" 2 , 
of which 13 rad m" 2 is due to errors associated with mea- 
surement and foreground removal, 12-17 rad m" 2 is due 
to residual small-scale fluctuations in the Galactic RM, and 
10-15 rad m" 2 is extragalactic Faraday rotation that is not in- 
trinsic to the emitting source and must originate in intervening 
systems between the radio sources and the Milky Way. 

We find a strong distinction between the fractional polariza- 
tions of radio sources whose optical counterparts are galaxies 
and those whose optical counterparts are quasars. The for- 
mer can be highly polarized, representing the extended lobes 
of radio galaxies, while the latter are only weakly polarized, 
representing the radio cores of active galaxies. Beyond this 
bimodality, we find no evolution of polarized fraction with 
redshift for either the galaxy or quasar population considered 
separately. 

Finally, we identify a strong depolarization effect in our 
RRM-redshift catalog, whereby sources with even modest 
residual Faraday rotation, |RRM| > 20 rad m~ 2 , have substan- 
tially reduced polarized fractions compared to sources with 
RRMs near rad m" 2 . We interpret this as beam depolariza- 
tion due to small-scale fluctuations in magnetic field strength 
and gas density in the same intervening population that con- 
tributes to the RRMs. 

A full consideration of the nature of the extragalactic source 
population that produces the observed Faraday rotation and 
depolarization is beyond the scope of this paper, but future 
investigations should consider the relationship of RRM with 
z for different sub-populations, and should study the depen- 
dence of the observed depolarization on wavelength, observ- 
ing frequency and angular resolution. Our analysis highlights 
the limitations imposed by the contribution of the Galactic 
foreground RM, even in cases where the foreground is mod- 
eled with a sophisticated algorithm that uses more than 40 000 
RMs as input. A much denser RM grid, with each RM deter- 
mined far more robustly than from the two frequency channels 
of NVSS data used here, is required to accurately account 
for the foreground contribution and to identify subtle trends 
of polarization properties with redshift. Such a data set will 
be provided by the upcoming Polarization Sky Survey of the 
Universe's Magnetism (POSS UM) on the Austra lian Square 
Kilometre Array Pathfinder dGaensler et al.ll2010l) . This next 
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generation of polarization data will be a powerful discrim- 
inant between different mechanisms for extragalactic Fara- 
day rotation measure and depolarization, and can thus pro- 
vide sensitive probes of magnetic field and electron density as 
a function of cosmic time. 
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Table 1 

Association classes between NVSS radio sources o f Taylor et al. 1 2009) and their 
optical counterparts. 



Class Description 



A Match using FIRST 

A(i) .... Close core match with FIRST 

A(ii) .... Core match with FIRST 

A(iii) . . . Single NVSS; closely matched double in FIRST 

A(iv) . . . Single NVSS; double in FIRST 

A(v) .... Single NVSS; closely matched triple in FIRST 

A(vi) . . . Single NVSS; triple in FIRST 

A(vii) . . Double-lobed in NVSS with FIRST 

B Core match in NVSS; either no FIRST detection or no FIRST data 

C Manual visual identification 

D Association made in NED 

E Association made in SIMBAD 



Table 2 

Prioritization of surveys and 
databases in selecting a final 
redshift, in cases where a radio 
source from Taylor et al. 1 2009) 
has matches in more than one 
optical redshift catalog. 



Priority 


Survey or Database 


1 


SDSS 


2 


6dFGS 


3 


2dFGRS 


4 


2QZ/6QZ 


5 


NED 


6 


SIMBAD 



Table 3 

Selected Columns from Our RM-Redshift Catalog of Polarized Radio Sources with Optical Counterparts 



NVSS ID 


SDSS 


2dFGRS 


6dFGS 


2QZ/6QZ 


NED 


SIMBAD 


Selected Value 


RM 




z Class 


z 


Class 




Class" 


z Class 




z 


Source 




(rad irT 2 ) 


(1) 


(2) (3) 


(4) 


(5) 


(6) 


(7) 


(8) (9) 


(10) 


(11) 


(12) 


(13) 


(14) 


J000010+305559 














1.801 


1.8 


NED 


1.801 


-37.9 ± 11.0 


J000030-112119 








0.10285 


A(iv) 








6dFGS 


0.10285 


+6.4 ±12.4 


J000132+145609 


0.39878 B 












0.39892 




SDSS 


0.39878 


-34.9 ±3.8 


J000153-302508 








1.30252 


B 




1.30252 




6dFGS 


1.30252 


-0.1 ±6.2 


J000154+020453 














0.402 




NED 


0.402 


-14.9 ±10.5 


J000255-265451 




0.0667 


B 


0.06666 


B 




0.06659 


0.0665 


6dFGS 


0.06666 


+0.8 ±7.2 


J000322-172711 














1.465 




NED 


1.465 


-33.2 ±2.4 


J000325-272631 




0.2501 


B 












2dFGRS 


0.2501 


+1.9±12.3 


J000327-154706 














0.508 


0.508 


NED 


0.508 


-2.8 ±4.0 


J000342-115149 








1.30999 


B 




1.30999 




6dFGS 


1.30999 


-3.8±13.3 



Note. — Table|3]is published in its entirety in the electronic edition of The Astrophysical Journal Supplement Series. The first 10 rows are shown here for guidance 
regarding its form and content. 

a The high proportion of 6dFGS and B class detections compared to the catalog as a whole is due to the low RAs of this sample (cf. Figure[T|) 
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Table 4 

The contents of our RM-redshift catalog, compiled from four optical surveys and two databases, broken down by 

association class. 



Class 


SDSS 


2dFGRS 


6dFGS 


2QZ/6QZ 


NED 


SIMBAD 


All 


Selected 


Selected, 




















\b\ > 20° 




1 135 


j ^ 


93 


27 






1200 


1 162 


1 160 


A(i) 


528 


6 


1 1 








JO 1 


539 


538 


ACiil 


12 


() 


() 


1 






13 


12 


12 


A(iii) 


225 


3 


5 


3 






236 


231 


231 


A(iv) 


83 


2 


2 


2 






89 


86 


86 


A(v) 


220 


4 


3 


5 






232 


225 


224 


A(vi) 


17 





2 









19 


19 


19 


A(vii) 


50 















50 


50 


50 


1! 


77 


43 


236 


35 






391 


375 


345 


C 


174 


9 


19 


3 






205 


195 


193 


D 










3186 




3186 


2023 


1742 


E 












1907 


1907 


248 


210 


Total 


1386 


67 


278 


65 


3186 


1907 


6889 


4003 


3650 


Selected 


1380 


44 


264 


44 


2023 


248 


4003 


N/A 


N/A 


Selected, \b\ > 20° 


1376 


44 


234 


44 


1742 


210 


3650 


N/A 


N/A 



Note. — The row labelled Total' gives the total number of associations made between Taylor et al. ( 2009) sources and the optical 
survey or database in question, while the rows labelled 'Selected' indicate the number of case s where this survey or database is the source 
of the "selected redshift" (see j j3.3l — for example, there are many cases for which the same Taylor et al. 1 2009) source has a redshift in 
both SDSS and NED, in which case the SDSS redshift is select ed over the NED on e. The column labelled 'All' lists the total number of 
times a given association class is assigned to a match between a Taylor et al. 1 2009) source and an optical counterpart, while the columns 
labelled 'Selected' indicate the number of selected redshifts in each association class — for example, most SIMBAD matches (class E) 
also have another match in class A, B or C from an optical survey. 



Table 5 

Selected Columns from Our RRM-Redshift Catalog of Polarized Radio Sources with Optical Counterparts 



NVSS ID 


e 


b 


NVSS RM 


Stokes / Flux 


Frac Pol 


Source 


Class 


Object 


Redshift 


GRM 


RRM 




n 


(°) 


(rad irT 2 ) 


(mJy) 


(%) 






Type 




(rad irT 2 ) 


(rad irT 2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


J000010+305559 


110.15 


-30.66 


-37.9 ±11.0 


88.2 ±2.7 


6.5 ±0.3 


NED 


D 


Quasar 


1.801 ±0.007 


-74.2 ±7.5 


+36.3 ±13.3 


J000030-112119 


83.29 


-70.20 


+6.4 ±12.4 


80.9 ±2.9 


6.8 ±0.4 


6dFGS 


A(iv) 




0.10285 ±... 


+1.9±4.1 


+4.5 ± 13.1 


J000132+145609 


105.37 


-46.23 


-34.9 ±3.8 


314.6±11.1 


5.5±0.1 


SDSS 


B 


Quasar 


0.39878±... 


-28.9 ±5.2 


-6.0 ±6.4 


J000153-302508 


13.14 


-78.66 


-0.1 ±6.2 


173.8 ±5.2 


6.4 ±0.2 


6dFGS 


B 




1.30252±... 


+15.4±4.1 


-15.5 ±7.4 


J000154+020453 


98.76 


-58.45 


-14.9 ±10.5 


300.7 ±9.7 


3.3 ±0.2 


NED 


D 


Galaxy 


0.402 ±... 


-8.3 ±5.4 


-6.6±11.8 


J000255-265451 


31.30 


-79.20 


+0.8 ±7.2 


94.4 ±3.2 


12.8 ±0.4 


6dFGS 


B 


0.06666±... 


+4.9 ±3.9 


-4.1 ±8.2 


J000322-172711 


71.53 


-75.28 


-33.2 ±2.4 


2414.8 ±72.4 


1.5 ±0.0 


NED 


D 


Quasar 


1.465 ±0.003 


+1.9±4.6 


-35.1 ±5.2 


J000325-272631 


28.49 


-79.33 


+1.9±12.3 


144.0 ±5.1 


5.2±0.3 


2dFGRS 


B 




0.2501 ±... 


+5.9 ±4.1 


^t.0±13.0 


J000327-154706 


76.02 


-74.10 


-2.8 ±4.0 


527.3 ±15.8 


3.2±0.1 


NED 


D 


Galaxy 


0.508 ±... 


-5.4 ±3.4 


+2.6 ±5.3 


J000342-1 15149 


84.35 


-71.07 


-3.8±13.3 


351.7± 12.2 


1.6±0.1 


6dFGS 


B 


1.30999 ±0.0068 


-2.4 ±3.3 


-1.4± 13.7 



Note. — Table [5] is published in its entirety in the electronic edition of The Aslrophysicai Journal Supplement Series. The first 10 rows are shown here for guidance regarding its form and 
content. 




Figure 1. The sky coverage for each of the surveys used in the creation of our RM vs. redshift catalog, shown using sparse sampling. 
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Figure 2. RMs of 37 543 extragalactic sources over the sky above a declination of —40° from the survey by Taylor et al. (2009). Red circles correspond to 
positive RMs while blue circles are negative. The size of the circle scales linearly with the magnitude of the RM. The sources are mapped onto a Galactic 
co-ordinate system. 
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Take an optical 
source 



Match with Taylor et 
al. (2009)-out to 3' 



✓ \ 



3+ matches - 
inspect visually (Class C) 



1 match found 
with NVSS 



1 



Look at offsets of 
match to optical 



Close NVSS match 
(<10") -Class B 



Offset <3"- 
Class A (i) 



Other criteria 
Class A (ii) 



Offset <3"- 
Class A (iii) 



Doubles criteria 
Class A (iv) 



Offset <3"- 
Class A (v) 



Doubles criteria 
Class A (vi) 




Offset criteria not 
met- single match 



Match optical with 
FIRST to 30" 



No match to 
FIRST 



One match 
to FIRST 



Two matches 
to FIRST 



Three matches 
to FIRST 



2 matches found 
with NVSS 



1 



Look at offsets of 
matches to optical 



Offset criteria met- 
candidate doubles 



Match optical with 
FIRST to 30" 




Complex - inspect 
visually (Class C) 



Offset <3"- 
Class A (vii) 



4+ matches - 
inspect visually (Class C) 



Figure 3. A simplified flowchart representation of the algorithm implemented in this paper to make associations between optical and radio sources in various 
surveys. Final association classes are shown shaded, with the same shading used for classes with similar admission criteria. 
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Figure 4. Examples of class A(i) associations (close core match with FIRST data) made with our algorithm. Optical data are shown in gray-scale with NVSS 
radio contours overlaid in red and FIRST radio contours overlaid in blue, (a) NVSS J140445-013021: class A(i) association for both SDSS and 2QZ/6QZ. 
(b) NVSS JO 10644- 103409: class A(i) association for both SDSS and 6dFGS. All contour levels are measured in mJy beam -1 . In this and subsequent related 
Figures, optical data are taken from SDSS where available, and otherwise from SuperCOSMOS. 




Figure 5. As for Figure|4] but showing examples of class A(ii) associations (core match with FIRST data), (a) NVSS J162530+270544: class A(ii) association 
with the SDSS. (b) NVSS J094735+583048: class A(ii) association with the SDSS. 
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Right Ascension (J2000) Right Ascension (J2000) 



Figure 6. As for Figure[4] but showing examples of class A(iii) associations (single NVSS source, closely matched double source seen with FIRST), (a) NVSS 
J01 1013-021954: class A(iii) association with the 6dFGS. (b) NVSS Jl 12055+173854: class A(iii) association with the SDSS. 
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Figure 7. As for Figure|4] but showing examples of class A(iv) associations (single NVSS source, double source seen with FIRST), (a) NVSS J1041 12+004550: 
class A(iv) association with the 2QZ/6QZ. (b) NVSS J230545-003608: class A(iv) association with the SDSS and 2dFGRS. 
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Figure 8. As for Figure [4] but showing examples of class A(v) associations (single NVSS source, closely matched triple source seen with FIRST), (a) NVSS 
J125003-013226 : class A(v) association with the 2dFGRS. (b) NVSS Jl 15233+493937: class A(v) association with the SDSS. 




Figure 9. As for Figure|4] but showing examples of class A(vi) associations (single NVSS source, triple source seen with FIRST), (a) NVSS J095406-065957: 
class A(vi) association with the 6dFGS. (b) NVSS J1201 15+180934: class A(vi) association with the SDSS; here, the faint optical source lies within the second 
FIRST component from the bottom. 
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Figure 10. As for Figure [4] but showing examples of class A(vii) associations (double-lobed radio source in NVSS, also detected in FIRST), (a) NVSS 
J213515-005255 and NVSS J21351 1-005233: class A(vii) association with the SDSS; the optical source lies in the middle of the FIRST contours, (b) NVSS 
J144558+122228 and NVSS 1144602+122258: class A(vii) association with the SDSS. 




Figure 11. As for Figure|4] but showing examples of class B associations (core match with NVSS but either no FIRST detection or no FIRST observations), (a) 
NVSS J000935-321636: class B association for both 6dFGS and 2dFGRS. (b) NVSS J002430-292848: class B association with the 2QZ/6QZ; in this case, the 
optical source is slightly offset from the radio center. 
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Figure 12. As for Figure[4] but showing examples of class C associations (complex sources requiring manual visual inspection), (a) NVSS J132644+031158: 
class C association with the SDSS; here the complexity is the presence of four FIRST components, (b) NVSS Jl 13305— 040047: class C association with the 
6dFGS; the complexity here is the extremely extended nature of the NVSS source, which FIRST does not detect. 



500 



400 



300 



200 



100 




500 



400 



300 



200 



100 



J 


1 (b) \ 





-10 10 

(RA Tay ior - RA optkal ) [arcsec] 



20 



-20 -10 10 

(Dec Taylor - Dec optical ) [arcsec] 



20 



Figure 13. Histograms showing the distribution of radio and optical positional offsets in our RM-redshift catalog in (a) Right Ascension and (b) Declination. The 
catalog includes a small number of associations with offsets outside the range shown (these sources, like all those beyond an offset of 7", have a double-lobed 
or complex radio morphology, leading to large separations between radio and optical source positions). 
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Figure 14. Map of the foreground Faraday sky in rad m 2 by OoDermann et al. 12012). 
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Figure 15. Comparison of (a) RM vs. z and (b) RRM vs. z, showing the impact of the GRM subtraction. 
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Figure 16. A comparison of our RRM-redshift catalog (bottom panel) with previously published data sets. A catalog of 268 objects discussed by Kronb erg et al.l 
1 2008) has not yet been made available. 
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Figure 17. RRM vs. redshift for (a) the whole RRM-redshift catalog of 3650 sources, (b) the subset of 1376 sources with selected redshifts from the SDSS, (c) 
the 516 galaxies with selected redshifts from the SDSS, and (d) the 860 quasars with selected redshifts from the SDSS. In each panel, individual data points are 
shown in black, while the mean RRM and ±lcr values of RRM on either side of this mean in independent adjacent bins of 50 points are shown by the solid red 
line and the dashed red lines, respectively. Note that there are 15 sources with RRM < -100 rad irr 2 and 10 sources with RRM > +100 rad itT 2 that are outside 
the range plotted. 
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Figure 18. As for Figure [T7l but showing fractional polarization vs. redshift. 
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Figure 19. As for Figure 1171 but showing RRM vs. fractional polarization. There are 15 sources with RRM < -100 rad nT 2 and 10 sources with RRM 
> +100 rad m~ 2 that are outside the range plotted. 



